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A realization of the approximation of the SU(3) leading representation with the 
microscopic Hamiltonian and a nucleon-nucleon interaction is presented in detail. An 
effective Hamiltonian reproducing results of calculations with some known potentials 
is constructed. It is shown that its structure is quite similar to that of the triaxial 
rotator, and the wave functions in the Elliott's scheme are linear combinations of 
Wigner's D-functions although they should be properly normalized. 

o 

1 Introduction 

For the theoretical studies of excitation spectra of valence nucleons in light and medium 
nuclei, Elliott [p]] proposed the basis states of irreducible representations (A,//) of the SU(3) 
group satisfying the Pauli exclusion principle. A special attention was paid to the most 
^ ■ symmetric irreducible ("leading") representations. In the states belonging to the leading 
representations (LR) the number of even nucleonic pairs is the largest, as they are character- 
ized by an even orbital momentum of the relative motion of the nucleons in the pairs, while 
the interaction between these nucleons are reproduced by the components V2S+i=3,2T+i=i 
and V13 of a central nucleon-nucleon (NN) potential. The number of odd nucleonic pairs 
is, on the contrary, the smallest, the orbital momentum is odd, and the acting components 
of a NN potential are V33 and V\\. While the even components are describing the strong 
attraction between the nucleons, the odd ones are responsible for the repulsion at short 
distances necessary to provide the saturation of nucleon forces and the known dependence 
of the volume of the nucleus on the number of nucleons. However, in the Elliott's scheme, 
the main argument for the importance of the LR was related not to the properties of the 
nuclear force, but rather to the properties of the operator of quadrupole (QQ) interaction 

-77QQ = G 2 - iL 2 , 

which was chosen to describe the residual interaction in the system. Here, 77 is a positive 
parameter, G2 is the second-order Casimir operator of the su(3) algebra, and L is the 
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squared orbital momentum. The eigenvalue #2 of the operator G2 is largest in the LR. As 
a result, firstly, the minimal eigenvalue — 77QQ is found in the LR, and secondly, the model 
reproduces the ordering of lowest levels of the principle rotational band. 

The Elliott's scheme was further developed by the appearance of the pseudo-SU(3) 
model0, ||, || and especially by the studies of the scissors mode0. In these studies, 
it became necessary to improve the phenomenological Hamiltonian and the wave functions 
of even-even nuclei (the ground states and the final states of the isovector Ml transitions) ||, 
0, ||. One more possible application of the SU(3) model is the neutron-rich nuclei like 
9 ' n Li and 10,11 Be, where the calculations may be made with some microscopic NN potential 
instead of just a phenomenological Hamiltonian. One has to find such an effective Hamil- 
tonian expressed in terms of the SU(3) group generators, which would generate the same 
spectrum as the microscopic Hamiltonian. One has also to find those 0(3)-scalar combina- 
tions of the SU(3) generators which, in principle, can enter such a Hamiltonian. This work 
deals with this kind of problems. 

We shall see that the Elliott's scheme and the triaxial rotator model share many common 
features. In particular, the basis states in both models can be expressed in terms of the 
Wigner's D-functions, but the normalization differs due to the fact that the density matrix in 
the Elliott's scheme corresponds to a mixed state, as opposed to a pure state in the case of the 
triaxial rotator. This conclusion is supported by the structure of the effective Hamiltonian in 
the Elliott's scheme. It has a form of a linear combination of scalar expressions made of the 
SU(3) group generators, and can be expressed through integer powers of the Hamiltonian of 
the triaxial rotator. The larger the number of quanta in the valence shell, the largest power 
of the rotator Hamiltonian enters the Elliott's effective Hamiltonian. 

We apply the LR approximation to the 10 Be nucleus as an example. The NN interaction 
is simulated by the known Volkov||] and Minnesota[rO| potentials in order to test whether 
they can reproduce the observed excitations of 10 Be. Usually, the Elliott's scheme is used 
with relatively simple phenomenological effective potentials. We therefore show in detail 
the procedure in the case of a microscopic Hamiltonian which, of course, is a generalization. 
We construct the scalar expressions from the SU(3) group generators and show that their 
linear combination can be reduced to the Hamiltonian of the triaxial rotator. Finally, we 
determine the parameters of the phenomenological effective potential making it equivalent 
to the microscopic one in the LR limit. 



2 Elliott's scheme and triaxial rotator 



Elliott's scheme (in its LR approximation) shares many common features with the theory 
of the rigid triaxial rotator pT]. It becomes clear in a space where the basis functions of 
both models are represented by the spherical Wigner's functions. The tool for the tran- 
sition to such a space is the construction of the function generating a complete LR basis 
in the form of a Slater determinant. As a result, the matrix elements of various operators 
are calculated easily. Operators involved are LQL and LQQL which have a simple alge- 
braic interpretation |T^] and consisting of the SU(3) group generators, as well as microscopic 
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operators of central, spin-orbit and tensor NN interactions. 

Consider the case of the 10 Be nucleus. The quantum numbers of basis states are, in fact, 
known, and we just check that our approach generates them correctly. At the same time, 
the explicit form of basis functions will be established, along with the proper normalization, 
which is crucial in microscopic calculations of spectra and transition probabilities. 

We first define the orbitals, distinguishing between proton and neutron configurations 
and restricting the basis with the minimal allowed number of oscillator quanta. Below we 
omit the spin-isospin quantum numbers for simplicity. 

Two proton (s— and p—) orbitals are 



.2\ v/2>^ / 



h^\'-i -777 exp I - T — I and lw (r) = -^expf-^-l, (1) 



^3/4 y 2 J ^ 17rv ' 7T 3 /4 y 2 / 

where the unit vector u is the first independent variable in the space where the Wigner's 
functions will be defined later. There are two protons with different spin projections in each 
of these two states. 

The remaining six neutrons are allocated in pairs in three states, 

Mr) = ^374 ex P ( -J J > Mr) = — v4~ ex P ( > ( 2 ) 



y/2 ([w X U|r) / r 



<>2u\ 



7T 3 / 4 (UU) 



{ - V ] • (3) 



There appeared another vector variable, the unit vector w, orthogonal to the vector u. 
We also introduce the vectors u and w for the conjugated orbitals. These vectors are also 
mutually orthogonal, but, in general, have a different orientation in space. 
It is easy to see that 



J M r )Mr)dr = 1, J 0^(r)0 lw (r)rfr = (uu) 
J (p* 0u (v)(p 0u (r)dr = 1, J #„(r)0i„(r)tfr = (uu) 



[r)dr 



WW 



uu 



Now we multiply these orbitals by the spin-isospin functions and construct the Slater 
determinant ^ of the nucleus 10 Be. A second Slater determinant \P* is constructed on the 
conjugated functions. It can be considered as a result of rotation of the coodinate frame 
transforming the vectors u, w into u, w. Therefore, the forthconing procedure is nothing else 
but a Peierls-Yocozz method of the angular momentum projection[^]. The overlap kernel 
(the result of integration of the product of two Slater determinants over all single-particle 
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vectors rj, % = 1. . .10) is expanded over the Wigner's spherical functions. The kernels of 
different operators are convenient to use instead of the wave functions since the number of 
variables is reduced. 

The overlap integral is actually very simple expression, 

J ***drxdr 2 • • • dr w = J( 10 Be) = (uu) 2 (ww) 2 . (4) 

The nucleus 10 Be has two protons in its p-shell, and the symmetry indices (A, /i) of the 
LR are A = 2, /i = 2. It follows the relation (|4]), which can be rewritten in another form due 
to the fact that the vectors involved are of unit length, 

/( 10 Be) = H 2 (ww) 2 = ^ 2 , (5) 

where du, 622 are the elements of the rotation matrix. These elements depend on the Euler 
angles only, so that (uu) = du, (ww) = 622- The Elliott's basis appears if the overlap 
(||) is expanded over the Wigner's /^-functions depending on the same Euler angles. This 
expansion is a necessary element if the Peierls-Yocozz method is implemented, because the 
weights of different angular momentum states have to be found. 
The overlap integral for a representation (A, n) is 

I(X,fi) = d x n d^ (6) 

Finding the coefficients C^. ^ of its expansion over the D-functions complete the analysis of 
this expression. First we write this expansion as follows, 

= E E C L KtR D L K k . (7) 

L K,K 

Evidently, g. = C^. K , so that the expansion matrix is Hermitian. 

One can see which D-functions enter the expansion ([?]) by invoking the concept of the 
point group D2 [M , elements of which alternate the sign of one of the vectors u, w (u, w), 
or the vector orthogonal to both of them. If the indices A and /1 are even, the expression 
(uu) A (ww)' 1 is invariant with respect to those transformations. Therefore it may contain 
only those D-functions which belong to the symmetric representation of D2, 

D° D 2 D 2 D 2 D 2 D 3 etc 

(ref. 0). 

Yet not all of these functions may enter the expansion. Its actual content depends on 
the choice of the coordinate frames used. The number of basis states does not depend on 
this choice, but the structure of the states does. We illustrate this on the example of the 
overlap integral (|5]). Following Elliott, we direct the z (z) axis along the vector u (u), and 
x (x) along w (w). Then the expansion takes the form 

d 2 d' 1 — D" I - I > _u_L_in- + -D 2 1 + 

u ll u 22 _ -| ^0,0 ^ o-\ \ A 0,0 ^ A -^0,2+ ^ ^ ±y 2+,0 ^ / ^ u 2+,2+ I ^ 
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0,2+ 



_^n4 _^n4 ,5 4 \ 

0,0 g -^0,2+ g 2+,0 g 2+,2+ I ' W 

The five coefficients at the .D-functions are the weight factors of the basis states, and their 
sum equals 1 as it should be. Moreover, the expansion (||) means that the Z)-functions 
are maps of the basis states ip1 m , where m is the projection of the orbital momentum on 
the external axis, a is an additional quantum number which may be necessary. Thus, for 
the ground state of 10 Be, i[)q — 1 with the weight 2/15, and for the state with L = 3, 
ipsm = D\_ m with the weight 1/6. One of the states with L = 2 has the weight 5/21 and 
the map 

ib x --D 2 +—D 2 
while the other, with the weight 1/3, 

,/,2 _ yj[ D 2 _ I n2 

r2m — 2 0,m 2 2 +' m ' 

Finally, the state with L = 4 has the weight 9/70 and the map 



2 o \/5 



m 



The wave functions obtained here are identical to those of the non-axial rotator model |T I 
the case when the non-axiality parameter 7 = tt/6. Thus we have established a link between 
these two models which are based on essentially different suggestions. 

The overlap integral, both in the general (|7|) and in the particular (||) cases, is also a 
density matrix calculated by integrating of the product of two Slater determinants over all 
single-particle variables. Unlike in the standard density matrix |0| where the integration 
is performed over some of the single-particle variables while the remaining ones become 
independent parameters, here the independent parameters are the vectors u, w (u,w). In 
other words, the introduction of the density matrix is accompanied by a transition from the 
space of single-particle variables to the space of Euler angles. In this new space the basis 
states of SU(3) representations are elegantly expressed as the Wigner's spherical functions. 

The density matrix is diagonal in the basis defined above. At the same time, the basis 
functions differ in their weight factors. That is why the density matrix describes a mixed 
state of the nuclear system, because otherwise the weight factors would be equal for all the 
functions. 

Let A and \i be even. Then the following expansion holds, 

i{\ ^) = E£ aO E r L J a Lm , (9) 
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where w%(\, /i) are the weight factors satisfying the condition 

£X>2(A,Ai) = l. 



If A > /i and L is even, then 



If L is odd, 



In any case, 



L a 



^Lm — ^2 ^LK^K+,m- 



K=0 



K=0 



(10) 



(11) 



(12) 



In particular, 



E^ia /ret' s: 
^LK^LK — °a,a'- 



K=0 



W (X,fl) 



(A- !)!!(// -1)!!(A + //)!! 
A!!/i!!(A + /i + l)!! 



(13) 



(14) 



2V ' w V ' w 2 (A + 2)( / u + 2)(A + /i + 3) \V3 A-// V J 



°20 



v / TTy ? + l 
\| 2v^+F : 



r 1 

°22 



\| 2v/T+F ' 



where 



= ^ (A + /i + 2)(A-AQ 
^ V A(A + 2)+ yU (/i + 2)' 



(16) 



5 A(A + 2)+ / x( / i + 2) [ 2 A/i + A + /i /— — ] 

2 (A + 2)( At + 2)(A + /x + 3) \V3 A-// y V ^ J ' 1 ' 



r 2 



\l 2v/T+y 



2 > ^22 



r 1 - - 



v^+?+ 1 

\| 2v^+F ' 



(18) 



The actual number of terms in @ is defined by the Elliott's rule for the basis states of an 
irreducible representation of the SU(3) group. 
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Elliott's choice of the axes simplifies the classification of the basis states. However, 
its drawback is a false impression that the Elliott's scheme contradicts the experimental 
evidence related to the structure of the wave functions. In the case of 10 Be, there are two 
D2-symmetric states with L = 2, and in neither of them the projection of the angular 
momentum on an internal axis K is a constant of motion, although experiments show that 
the rotational states are grouped in bands K — 0, K = 2, etc. 

In order to solve this contradiction, we redirect the axis z (z) along the vector orthogonal 
to u and w (u and w). This results in a different structure of the expansion with the same 
weight factors, 

d 2 d 2 - — D° + — D 2 + - D 2 + - D 3 + 

U 11 U 22 — "I r "^0,0 ^ 91 "^0,0 ^ n- U 2+,2+ ^ n 1J 2-,2-^ 



+ ~ { ~D\ +A+ - y^Dj +0 - ^DU + + ~D° 00 ) • (19) 
70 \36 4+ ' 4+ 36 4+ '° 36 °' 4+ 36 °'°J V ; 

One can see here that this choice of the rotation axis conserves K if L — 2 which makes a 
better agreement between the Eliott's scheme and the experimental data. As for the state 
with L — 4, its main component K — 4 enters with the weight factor 35/36, the component 
K = 2 is missing, and the component K = 2 has a small weight of 1/36. Such a drastic 
difference in the structure of states L = 2 and L = 4 may be related to an experimentally 
observed phenomenon; in the first rotational band, the value of K undergoes a sudden 
change when L reaches a critical value. 

In general, the best choice of the rotation axis may be made according to the following 
rules, depending on the SU(3) symmetry indices (A,//) of the LR. 

1. If the indices of the LR are (A, 0), the overlap integral is 

(uu) A , 

and the rotation axis should be directed along the vector u (u). Then, all allowed 
states have K — 0. If /i A, components with K 7^ appear, but their amplitude 
is small in the states of the main rotational band, so that they can be approximately 
treated as K = states. 

2. Similarly, if the indices of the LR are (0,/i), the overlap integral is 

(wwf , 

and the rotation axis should be directed along w (w). Again, all the states have K — 
only, and if A fi, K may be considered approximately equal to zero. 

3. In the X — pL case, the rotation axis should be directed along [u x w] ([u x w]). If 
A and ji are slightly different, a small admixture of K = 2 states appear in the main 
rotational band even at L = 2, and K — admixtures are found in the K = 2 band. 
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3 Overlap integral with the Hamiltonian 

The next stage in the realization of the Elliott's scheme for the 10 Be spectrum is the cal- 
culation of the overlap integral of the generating Slater determinant with the Hamiltonian 
H = T + U, where U is a central exchange potential having a Gaussian form, 

= A{h + I 2 ) + Bh + (P~ 2A)I, (20) 

where 

I t = M V 2 (ww*) 2 , I 2 = wV 2 (uu*) 2 , h = ([u x u*][w x w*])(uu*)(ww*). 

Below we leave only those terms which split the spectrum and, therefore, are of interest. 
The coefficients A, B and P are expressed through amplitudes of the even components V31 
and V13 of the nucleon-nucleon potential, and 

1 2r 2 

K 

where r is the oscillator length, b is the range of the potential. 

A = (^f) V 13 , B = z^ (l^) 2 (Via + 3V 3 i). (21) 

If we direct the axis z (z) along [u x w], ([u x w]), we obtain 

2 1 2 2 

/1 + J 2 = g A),o + 3-^0,0 + D 2+,2+-, 





0,0 



Js ~ 30^°°° + 21 jD ° 2 ° 70 l36^ + ' 4+ 36 ^ 36 D *>*+ + 36^' 



It is easy now to find the energies of the five states of 10 Be. The ground + state is 
located at 

E = 5A + ^B. (22) 
E will be used as a reference point for the excited states, 

Bt -£ = -^(A-f)- \b, (23) 



El - E = -2 [ A - 



B\ 3 



4/ 4 



E3 — E Q — 



-5U-f 



B, 



5* 



(24) 
(25) 



E, MeV 



6.9 



L = 4- 



E, MeV 
- 7.25 




Figure 1: Spectrum of 10 Be for the Hamiltonian ( p0|) with Volkov and Minnesota interac- 
tions. 

Ek-Eo = -s(A-j^-^B. (26) 

Note that if the Volkov potential is used, the states L = 2, K = and L = 2, K = 2 are 
degenerate, simply because this potential has V31 = V13. Fig. 1 shows the spectrum of 10 Be 
for the Volkov potential with m = 0.6, if ro = 1.64 fm. This value of ro is optimized by 
variation; it minimizes the energy of the ground state and correspond to the r.m.s. radius 
2.29 fm (the experimental value is 2.3 ± 0.2 fm[|i~6|l). 

While the Volkov potential is consistent with the observed value for the radius of 10 Be, 
the Minnesota potential with u = 0.98 (this value of the exchange mixture parameter seems 
to be the best choice for the nuclei in the first half of the p-shell) leads to ro=1.43 fm, 
and besides, yields a wrong sequence of excited states. The saturation is not reached for 
both potentials but, again, this is more evident in the Minnesota case (light nuclei being 
considered). Nevertheless, if r is set to 1.64 fm in the Minnesota case, one obtains the 
spectrum shown in Fig. 1 on the right. Now the 2 + states split as V31 7^ V13. However, the 
splitting is just 1.29 MeV, which is twice smaller than observed (2.59 MeV). The splitting 
E\ — E\ is proportional to V31 — V13. Therefore, this value has to be increased in order to 
achieve a good description of the experimental results within the model. This suggestion is 
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also supported by the fact that the value 



-{El + E\) -E = 1.125 MeV 

appears to be four times less than the experimental one, and almost twice less than in the 
Volkov potential case (2.07 MeV). 

Historically, both potentials were designed to describe the experimental data for the 
nuclei in the beginning of the p-shell. It is therefore not surprising that the their application 
to heavier nuclei may be accompanied by some corrections in the potentials themselves. 

The Elliott's scheme predicts the existence of the levels 3 + and 4 + . Energy of the former 
is close to the sum of the energies of the 2 + states. The wave function of the latter has 
predominantly K = 4 which affects the probabilities of E2 transitions. The isoscalar electric 
transition from the 4 + state to the L = 2, K = state is forbidden, leaving allowed only 
the transition to the L = 2, K = 2 state. 

It is important that a right choice of the central nucleon-nucleon interaction makes the 
states L = 2, K = and L = 2, K = 2 uncoupled. 



4 Effective Hamiltonian of the SU(3) model 



When the overlap integral ( |2"0"|) was calculated, the operator H was acting on the generating 
Slater determinant Its SU(3) symmetry was (2,2), and the result was projected to the 
states with the same symmetry. The projection was performed by an integration. Therefore, 
the overlap integral (^) keeps us within the basis of the representation (2,2). In order 
to understand how the Hamiltonian H acts in this representation, we define an operator 
H(u, w) as follows, 

(V\H\V) = H(u,w)(V\V). (27) 

All operators which acting onto the basis states of an irreducible representations produce 
only those states which belong to the same representation are known: they are the group 
generators. We are interested here only in those expressions constructed from the SU(3) 
generators which are scalars with respect to the 0(3) group. They are known, too [jP7|1 , 



lq| . They are the operator QQ introduced by Elliott, and the operators = LQL and 



Qi = LQQL fl2" |. All other operators may be written as polynomials of these three. 

In our model the wave functions are represented by superpositions of the D-functions. 
Therefore, as in the case of the Hamiltonian of the triaxial rotator, all operators may be 
expressed via the SU(3) symmetry indices A,/i and the three projections of the orbital 
momentum L to the internal axes: Li,L 2 ,L 3 . It has been shown 111 91] that 



- Q = aL{ + bL 2 2 + cL\ + fi 2 , (28) 



Six = a 2 L\ + h 2 L\ + c 2 L\ 



2 -V-Ll+ 1 -(LlLl + LlLl) + 2bQ 2 , 



(29) 
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where 



2A + /i + 3 , -X + n -A-2/X-3 , on , 
a = , o = , c = . (30) 



f^2 — g \LlL2L3 + L2L3L1 + L3L1L2 + + L3L2L1 + L1L3Z/2J . (31) 

The structure of the operators f2 and f2i is similar to that of the Hamiltonian of the 
triaxial rotator, but the inertial parameters (the coefficients of L 2 , L 2 , £3) are not the 
principle values of the tensor of inertia, but rather the principle values of the tensor of the 
intrinsic quadrupole momentum (in the case of ti) or their squared values (in the case of 
tii). Additional terms 

ti 2 and ^(L\Ll + LlLl) 

only lead to the finite number of basis functions of an SU(3) representation whereas in the 
triaxial rotator case the basis is infinite. 

The operators ti and tii commute with L 2 , therefore the orbital momentum L, as ex- 
pected, is a constant of motion. Besides, these operators are D2-invariant, so that their 
eigenfunctions have a definite D2 symmetry. However, among all eigenfunctions of ti and 
Cli, only those are of physical interest which have a D 2 symmetry identical to that of the 
basis states produced by the overlap integral (0). 

It is natural to seek the eigenfunctions <pL,m of ti (or ti\, or a linear combination of ti 
and tii) in the form of a superposition of the functions 

taking into account their normalization, that is 

H^ = Y. B ^4^tm- (32) 

a 

Then the coefficients Bl >0 , must satisfy the following set of equations, 

a\ti\L, a') - 5 a>a/ )B L)a , = 0. (33) 

a' 

One may calculate the matrix elements (L, a\ti\L, a'), by calculating first the matrix element 
of ti between the generating functions and then expanding the result, 

= EE y/wt(X,fi)wt(X lfl )(L, a\ti\L, a') £ ^J^ m . (34) 

L a,a' m 

Alternatively, a simple formula ( p5| ) for ti provides for another way, 

ti^lplm = £<L, "M ^ «')\/<<- (35) 
a' 
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The action of the operators Li, 1/2,1/3 on the D-functions is known (cf. ||15|| ). Both ap- 
proaches lead to the same Hermitian matrix, 

\\(L,a\n\L,a'}\\. 

Diagonalization of this matrix yields the eigenvalues of Q. At a given L, the size of the 
matrix depends on the number of possible values of a. 

Finally, one more possibility is to expand <pL,m over the .D-functions straightforwardly, 



N L £ A L , K D L K±>m . (36) 



K 



This way is the simplest, it results in a non-Hermitian matrix 

\\(L,K±\Cl\L,K'±)\\, (37) 

because (a term in both Q and Cli) contains an imaginary unit i as a factor, and that 
makes it a non-self-conjugate operator. This does not mean, however, that the matrix will 
necessarily have complex eigenvalues Q. The coefficients A LjK are found from the set 

]T((L, K ± \Q\L, K'±) - 8 K , K ,)A L)K = 0, (38) 

K' 

where the matrix elements (L, K ± \ fl\L, K'±) are defined by 

&D L K±tm = Y,(L,K± |0|L, K'±)D L K , ±>m . (39) 

K' 

The normalization factor Nl of the eigenvalues ([36]) is found from the condition 

^E4iEV^V = i. (40) 



• J W L K 

The coefficients C£ K were defined above ( (|TT| ) and ( |T2] ) ) . They can be found by expanding 
the overlap integral (^) over the basis states of an irreducible representation of the SU(3) 
group. The summation over K in (|D) signifies the projection of an eigenfunction of Cl to 
the allowed states. Hence those eigenfunctions of Cl for which each of the sums 

k 

vanishes are forbidden. They do not satisfy the Pauli exclusion principle and must be 
disregarded. With L increasing, such states appear, sooner or later, and this must be taken 
into account in the calculations of spectra. 



^ Hermitian matrix has real eigenvalues only, but a reverse statement is not right, in general. 
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5 General remarks on the effective Hamiltonian 



At least at the first stage, the effective Hamiltonian is better to be constructed as a linear 
combination of the operators QQ, Cl and Cli. It is now clear that the result will be, in 
fact, the Hamiltonian of the triaxial rotator model with some additional terms which may 
be important in the actual calculations of the LR spectra. The coefficients of the three 
operators may then serve as phenomenological parameters. 

The question is, can such a Hamiltonian describe the observed low-energy spectra of 
even-even nuclei? One may believe that the phenomenological parameters would be the 
same for a whole range of nuclei while the SU(3) symmetry indices would change from 
nucleus to nucleus according to known general rules. 



Such an approach has been used in |T7j and Ul8| , and there an attempt was made to 



extend the model to excited states with higher L. Then, one may question, whether it suffices 
to have the three operators entering the Hamiltonian in their first order, or a higher-order 
terms are necessary? 

In Ref. [EDI the authors studied the 20 Ne and 44 Ti nuclei and concluded, that the effective 



Hamiltonian of their LR (A, 0) ((8,0) for 20 Ne and (12,0) for 44 Ti), derived from the overlap 
integrals of the generating Slater determinants with the Gaussian interaction, contains first 
two orders of QQ for 20 Ne, and first three orders for 44 Ti. It became clear that higher orders 
are necessary for higher oscillator shells. The influence of higher-order terms is increasing 
with L. For representations (A,0), the operators Cl and f2i degenerate and reduce to the 
operator L 2 . Hence to derive the terms containing high orders of Cl and Cli, one has to use 
the LRs with A 7^ 0, fi 7^ and go beyond the p-shell. Alternatively, the higher-order terms 
may be included in the Hamiltonian phenomenologically. 



6 The mass quadrupole operator 

We now show that a, b, c are the principle values of the quadrupole momentum. To 
start with, the explicit form of the operator of the mass quadrupole momentum Q a @ in the 
intrinsic coordinate frame is shown in The same operator in the laboratory frame (with 
our choice of axes) is 

l ' •> • \ , V^/x , „\n2 ^* / n-> / ;n2 ; r>2 



Q 2m = 3 (-A + //)£H >0 + ^-(A + fx)Di >2+ - -y (L>l il+ Li + + D^ 2 _L 2 ) . (41) 



Here 



1 V3 
-(-A + //) and — (A + /x) 



are the principle values of the traceless tensor of the quadrupole momentum. 
It follows from (|l]) that 

QQ = ^(A 2 + A / u + / u 2 + 3A + 3 / u)-^L 2 . (42) 
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Also note that 

= iplnA+^m 1 A+ ~ ^m,X-^m' ,X- + ^m,2-^m',2-) • (43) 

We have arrived to an expected, D 2 -invariant expression. Indeed, acting on such an expres- 
sion, the operator of the quadrupole momentum must yield another D 2 -invariant expression. 
Consider now the Eliiott's choice of axes. 

Q2m = \ l (a - l -^)D 2 m>Q + ^-(b - c)D 2 m>2+ ^ + ... = i(2A + ^)D 2 mfi + ^D 2 m>2+ + ... 

By definition, 

a = ~(2A + /x) = a-l, b = ~(-A + n) = b, c = ~(-A - 2/x) = c + 1. 

o o o 

Then, for the representation (A, 0) one obtains 

2A 

Q 2m = y< + - (44) 

In this case, the axially symmetric nucleus has a prolate shape and therefore its intrinsic 
quadrupole momentum^ is directed along the rotation axis and is positive. 
Yet another choice, 

Q 2m = | {(c - ^)<o + %a - 6)< 2+ } + .- = 

= ~(2^A)< 0+ |a< 2+ + ... 
is convenient for the limiting case of (0, /x), when 

g 2m = -^< + - (45) 

Now the nucleus has an oblate shape, the quadrupole momentum is again oriented along 
the rotation axis but this time, negative. 

Finally, we show the matrix elements of the isoscalar El transition from the ground state 
to the states ip 2m °f the Elliott's basis when the indices A,/i are both even. 



^{^(-A + ^ 2 1 o + ^(A + ^ 2 1 2 } 



( L = 2 , a = 1, m\Q 2m \L = 0) = { -(-A + ^C*, + ^-(A + /i)^ ^ , (46) 
(L = 2, a = 2, m|Q 2m |L = 0) = ^(-A + //)C* + ^(A + /i)Cf 2 ) . (47) 




2 The quadrupole momentum is a traceless tensor with two main components, Q20 an d *?22+- In the 
axially-symmetric case the second component vanishes, and the quadrupole momentum is said to be oriented 
along one of its principal axes. 
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7 Effective Hamiltonian for Be 



We shall now try to built an effective Hamiltonian for 10 Be from L 2 , Q, &i. The only 
condition we impose on the Hamiltonian is the equivalence of its spectra (eigenvalues and 
eigenvectors) to that found above in the LR approximation. 

H(u, w) = B (J - i L 2 ) + [A - |) (p + qL 2 + rfii + sO) , (48) 

where p, q, r, s are the coefficients to be found. 

We have noted earlier that the matrix elements of the norm and Hamiltonian operators 
between the basis functions ip\ m and ip^m with L = 2 are diagonal only. Meanwhile, the 
eigenfunctions of Cl are linear combinations of these states. Therefore, s = 0, and it is the 
operator Cli that splits the two levels with L = 2. It follows (|23|) and (p4|), that 

El -El = \ (.4 - |) . (49) 

The eigenvalues of &i in these states are 

Q^L = -44 ^L, ^L = -12V4r (50) 

Using 

E = 5A + ~B, and £ 2 X -E = ~(A-^j- ^B, (51) 
we define p, q, r and obtain the effective Hamiltonian in its final form, 

AKw) = B (§-iL») + (^-f)(5-lf + lft l ). (52) 

One may check straightforwardly that this Hamiltonian yields the same values of E 3 and 
E4, too. 

In the phenomenological approach, the coefficients A and B should be chosen so that to 
reproduce the experimental spectrum of 10 Be. It would suffice then to use the energy of the 
first three states, £"0, E\ and E\. The energies E 3 and E^ could serve as a test. But their 
experimental values are not known at present. 

8 Summary 

We have shown that the realization of the SU(3) LR approximation with the microscopic 
Hamiltonian and a nucleon-nucleon potential (such as the Volkov and Minnesota potentials) 
may be reduced to the calculations with an effective potential. The latter is a linear com- 
bination of the operators L , Cl and Cli, as well as their higher orders. In many ways this 
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Hamiltonian is similar to that of the triaxial rotator model. However, the inertial param- 
eters of such a rotator are not inversely proportional to the principal values of the tensor 
of inertia, but proportional to the principal values of the tensor of the quadrupole momen- 
tum. Both the basis functions of an irreducible representation of the SU(3) group and the 
basis functions of the rotator model are superpositions of Wigner's D-functions, but their 
normalization should be calculated separately. 

Some problems appearing in the calculations of the spectra of nuclei in the second half 
of the p-shell with the Volkov and Minnesota potentials are discussed. 

Some remarks on the structure of effective Hamiltonians of medium and heavy nuclei 
are also given. 
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